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We examine the class of barotropic fluid models of dark energy, in which the pressure is an 
explicit function of the density, p = f{p). Through general physical considerations we constrain 
the asymptotic past and future behaviors and show that this class is equivalent to the sum of a 
cosmological constant and a decelerating perfect fluid, or "aether", with wae > 0. Barotropic 
models give substantially disjoint predictions from quintessence, except in the limit of ACDM. They 
are also interesting in that they simultaneously can ameliorate the coincidence problem and yet 
"predict" a value of w ~ —1. 



INTRODUCTION 



Observational evidence strongly points to an acceler- 
ated expansion of the universe [l|, but the physical 
origin of this acceleration is unknown. Since general rel- 
ativity relates spacetime curvature (and hence accelera- 
tion) to energy, it is natural to hypothesize that either 
this relationship must be modified (as in extended grav- 
ity models [1, or that there is some additional source 
of energy density driving the expansion. Many models in 
the latter category have been proposed, including scalar 
field models ("quintessence") [1, i, 0, H, i, 13 > 

scalar 

field models with modified kinetic terms ( "k-essence" ) 
m [H, [11 Q m, [M [13, [S [S HI, and even more 
exotic possibilities (see, for example, the review of (HI). 

Here we examine in more detail a class of models called 
barotropic fluids, in which the dark energy pressure pde 
is given as an explicit function of the density pde- This 
can be viewed as a very simple prescription, in contrast to 
other models (such as the ones above) where the relation 
is implicit, written in terms of intermediate variables. Al- 
though some barotropic models are well studied, we find 
here a number of new, general properties that makes this 
class of interest in dark energy physics. Specific models 
investigated previously include the Chaplygin gas [H, 
and the generalized Cha plygin gas [22, [2J|, the linear 
equation of state [1^, [H, \2T\. |28|| and the afRne equation 
of state m, [lOl (note these are actually the same model) , 
the quadratic equation of state [1^, [3r| , and the van der 
Waals equation of state [H, [s^ . Such models have been 
considered either as unified models for dark matter and 
dark energy together, or as models for the dark energy 
alone. We confine our attention to the latter case. 

One of the key properties of barotropic fluids is that 
the sound speed, = dpoE/dpoE, does not have to 
equal the speed of light as in quintessence models. In 
addition, the condition > causes the barotropic dy- 
namical behavior to be distinguishable from quintessence 
- that is, they tend to lie in distinct regions of the equa- 
tion of state phase space (sj). Here we extend these re- 
sults to a more general discussion of the types of behav- 



ior that are allowed for barotropic models. We will see 
that limits < < 1 on the sound speed, along with 
some fairly general observational constraints, allow sur- 
prisingly broad conclusions to be made about the proper- 
ties of viable barotropic models, ruling out some models 
in the literature. Moreover, the properties describe an at- 
tractively simple picture of dark energy, together with a 
possible simultaneous resolution of the coincidence prob- 
lem and why today the equation of state w is near —1. 

We discuss general properties of barotropic fluids in 
i flll together with some special cases. In i jllll we explore 
the distinction between barotropic and scalar field so- 
lutions to the dark energy puzzle, and the relation to 
ACDM, finishing with a comparison of the coincidence 
and w —I behaviors of the different classes. 



II. PROPERTIES OF A BAROTROPIC FLUID 

A. General Properties 

We define a barotropic fiuid as any fiuid in which the 
physics of the fiuid is fully determined by the pressure as 
an explicit function of the density: 



Pde = Upde)- 



(1) 



Thus, the equation of state function / completely char- 
acterizes a barotropic fluid. For example, the generalized 
Chaplygin gas has the equation of state function [23 | 



Pde 



A 

a ' 

Pde 



(2) 



where A and a are constants, while the quadratic equa- 
tion of state is 



Pde = Po + apDE + PPde^ 



(3) 



with the linear (or affine) model corresponding to the 
special case /3 = 0. The van der Waals equation of state 
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is |32|.|35 



IPDE 2 



(4) 



How do barotropic models differ from quintessence 
models for dark energy? For a canonical, minimally cou- 
pled scalar field 0, the relation between pressure and den- 
sity is given parametrically by 



P0 = + v 



(5) 
(6) 



where 4) = d(j)/dt, and V{4)) is the quintessence potential. 

One can certainly find quintessence potentials for 
which p0 is not a single-valued function of and which 
therefore can never be described in terms of equation llj . 
For instance, a field oscillating about the minimum of a 
potential has zero kinetic term at either extreme of the 
oscillation and zero potential term at the minimum (in 
terms of the equation of state ratio w — p^/pcf,, this is 
w = —1 and respectively). Thus, passes through 
twice on every oscillation, with a decreasing value of 
each time (since diminishes as the universe expands). 
On the other hand, many quintessence models can be 
characterized by a pressure which is a single- valued func- 
tion of density, so what distinguishes these models from 
the barotropic models we consider here?^ 

The main distinction is that barotropic models have a 
value of dpDE / dpDE which is constrained by limits on 
the sound speed. The sound speed for a barotropic fluid 
is given by 



dp 



DE 



dpDE 



(7) 



Note that there is no need to write partial derivatives, 
since pde depends only on puE- To ensure stability, we 
must have > 0, so that dp^E I dpoE > 0. Thus, the 
function /(pde) in equation ([T]) is not arbitrary; it must 
satisfy df /dpoE > 0. We will further require, for causal- 
ity, that Cg < 1 (d^, but also see [1^). This imposes the 
additional constraint df/dpjjE < 1- 

In contrast, a canonical minimally coupled scalar field 
is an imperfect fluid. While its adiabatic sound speed 
is given by equation ([7]) (written as a partial derivative 
holding the entropy, or scale factor, fixed), its physi- 
cal sound speed is always equal to the speed of light. 



^ Indeed, it is trivial to write down an effective potential for a 
barotropic model: 

= (p-f)/2 



/' 



and from this V{(t>) alone one could not tell if the physics was 
barotropic or quintessential. 



Thus, these models are not subject to the constraint 
that dp/ dp > 0; in fact, they generically have dp/ dp < 
[13, [33 • This is the reason that barotropic models and 
many quintessence models occupy disjoint regions in the 
w — w' phase plane [34i |. (For more general discussions 
of the behavior of scalar field dark energy models in the 
w-w' phase space, see [H [H IH, E^l ) • 

Starting from the definition of the equation of state 
ratio, w — pde/pde, and taking the derivative with re- 
spect to the logarithmic scale factor In a, denoted by a 
prime, one has 



w 



'3(1 + w) 



dp 



DE 



dpDE 



= ~3{1 + w){cl ~ w) . 
The requirement that > then gives [li 
w' < 3w(l -I- w). 



(8) 
(9) 

(10) 



Here we consider only barotropic models for which w > 
— 1, although it is also possible to gen erate barotropic 
phantom models with w < —1 (ill. |42|. [isj . Since, for 
dark energy, w < 0, we have w' < for all barotropic 
fluids that can serve as dark energy. Models of this kind, 
in which w approaches —1 with time, have been dubbed 
"freezing" models [s?! , although for quintessence freezing 
models one frequently has the opposite of equation (|10p : 
w' > 3w{l + w). (Quintessence models can be found 
to violate this, while barotropic models will never break 
Eq. [To]). Further, the upper bound on the sound speed. 



(11) 



< 1, gives a lower bound on w': 

w' > -3(1 + w)(l - w). 



This is precisely the null line for the w ~ w' phase plane 
[ssl ] and leads to the one exception where quintessence 
models are exactly equivalent to barotropic models: skat- 
ing models [H, |3| , with kinetic energy along a flat po- 
tential, follow the equality in equation ()11|1 and corre- 
spond to barotropic models with /{pde) — Pde — P*', 
both have — 1. This is a pathological case, however, 
as the kinetic energy of skating models redshifts as 
and the model quickly becomes indistinguishable from a 
cosmological constant. 

Since w decreases with time, equation ([H]) has a generic 
future attractor at w = —1, independent of the functional 
form of /{pde)- At the attractor, the density, pressure, 
and sound speed asymptotically approach constant val- 
ues, which we denote poo , Poo , and Cgoo • Taking Cg = Csoo 
in equation ([9]), we see that w approaches —1 as 



1 



(12) 



and the dark energy density asymptotically approaches 



Poo as 



Pde — Poo 



-3(l+c^^) 



(13) 
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Since at late times the density always approaches a 
constant, Poo, this suggests decomposing the barotropic 
energy density into two components, 

PDE = Poo + PAE , (14) 

where the first term represents an always present cosmo- 
logical constant and the second term defines an "aether" 
density, pae, a spatially pervasive fluid with sound speed 
generally different from unity. 

The aether component is itself barotropic (with zero 
cosmological constant piece), and by comparing equa- 
tions (HH) and (fT51) to the usual behavior of some com- 
ponent X, px ^ a-3(i+u.^)^ Qj^g gggg ^Yisx < Wae < 1- 
That is, we can represent any barotropic fluid overall 
as consisting of a cosmological constant plus a positive 
equation of state perfect fluid. The sound speed of the 
aether is equal to the sound speed of the full barotropic 
component, a result that is easily proven using the for- 
mulas for summing components, equations (26) and (27) 
of [40], together with equation 

Note that the cosmic acceleration from the dark energy 
arises purely from the cosmological constant piece. In 
contrast to the barotropic dark energy as a whole, the 
aether component acts to decelerate and has a density 
that decreases at least as fast as the matter density, pM ^ 
, and no more rapidly than a stiff fluid, p ^ . 

As a specific example, consider the Chaplygin gas 
model, which is given by equation ^ with a — \. The 
full Chaplygin gas density evolves as 

PDE = ^JpIo + {p\,Efl- Plo)a^^, (15) 

with the future attractor at poo = v^, so the aether 
density is given by 



PAE ^ ^ pIo + {pj^Efi- Plo)a ^ - Poo ■ (16) 

The aether density decays as a^^ at early times (i.e. 
WAE — 0), while in the limit where the constant-density 
attractor is approached, pae decays as (i.e. wae = 
+!)■ 

Requiring a matter dominated era at high redshift con- 
strains the behavior of the aether component at early 
times. Since wae > 0, the energy density of this compo- 
nent will tend to dominate in the past. However such an 
aether component with a density greater than the matter 
density at high redshift violates a host of observational 
constraints [45]. Thus, we must have wae ^ as a — > 0. 
Furthermore, applying equation Q to the aether com- 
ponent, we see that as a ^ we have < wae so in 
this limit ^ 0. 

These arguments then allow us to put very general 
constraints on the behavior of the overall barotropic dark 
energy: it will always behave like a pressureless dust com- 
ponent at early times, and like a cosmological constant 
at late times. (Note that "late times" could mean the far 
future, rather than the present). In order to produce 



dust-like behavior at early times, the functional form 
of /{pde) must satisfy the constraint df/dpoE — > as 
PDE — > oo, as well as the previously-discussed sound- 
speed limits: < df/dpoE < 1 at all times. We discuss 
the observational implications of this further in ijllll 

For the barotropic models we have mentioned, these 
limits impose severe constraints. For example, for the 
quadratic and afhne models of equation ([3]), our limits 
impose f3 = and a = 0, so all acceptable dark energy 
models of this type reduce to the simple case p^E — con- 
stant. This corresponds to the special case where pae 
behaves exactly like dust at all times (see the next sec- 
tion). Similarly, the van der Waals model (equation |4]) is 
found to have unphysical behavior. For the generalized 
Chaplygin gas (equation [2]) , our constraints give a > 0, 
ruling out several extensions of the generalized Chaplygin 
gas |46j . 

The aether decomposition also provides a simple recipe 
for producing acceptable barotropic models. Specifically, 
every barotropic equation of state can be written in the 
form 

PDE ^ -Poo + gipDE - Poo), (17) 
where the function g is subject to the constraints 

0<dg{p)/dp<l, (18) 

5(0) = 0, (19) 
and, in the limit where p oo, 

dg{p)/dp ^ 0. (20) 

B. Special Cases 

Now consider some special cases of interest. When 
= 0, we have the previously-mentioned constant pres- 
sure model, characterized by p = —poo- The density in 
this case evolves as 

p = Poo + Ca~\ (21) 

i.e., this looks just like the ACDM model (though with 
an additional matter contribution) . This model has been 
previously discussed elsewhere. For instance, it corre- 
sponds to the a = limit of the generalized Chaplygin 
gas [l^l , and it is a special case of the "mocker models" 
discussed in [sll , in which it was noted that such models 
are characterized by w' = 3w(l -\-w), defined there as the 
constant pressure line. 

For constant but nonzero, we have 

p = Poo + Ca-3(i+c^) (22) 

These models are all observationally excluded as noted 
earlier, since they asymptotically dominate the expansion 
at high redshift. In particular, pure skating models have 
= 1 and therefore correspond to dg/dp = 1 at all times; 
such models then violate the condition in equation (j20p 
above. 
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C. Relation to k-essence 

Note that there is a one-to-one correspondence be- 
tween barotropic fluids and the subset of k-essence mod- 
els with constant potential, the so-called "purely kinetic" 
k-essence models. The latter models were first investi- 
gated in the context of inflation fll] and later proposed 
as unified models of dark matter and dark energy p^ . 
These models are characterized by a Lagrangian of the 
form: 



where X is 



P = F{X), 



(23) 



(24) 



and is the k-essence scalar field. The pressure in these 
models is simply given by equation (|23p . while the energy 
density is 



p = 2X{dF/dX) ~ F. 



(25) 



To convert between a kinetic k-essence model and a 
barotropic model, one inverts the function F to find 
X = F~^{p) and substitutes this into equation ([^ 
to convert p{X) into p{p). Inverting this delivers the 
barotropic equation ([T]). If one starts with a barotropic 
model, one inverts / to give p = f~^(jp) and interprets 
equation ((25l) as a differential equation to solve for F{X). 
To wit (cf. [£3) 



dF 



r\F) + F 



= MCX^/^) 



(26) 



where C is a constant. This one-to-one correspondence 
(assuming the functions are invertible) means that the 
trajectories for purely kinetic k-essence lie in the same 
region of the w — w' phase plane as the trajectories for 
barotropic fluids [20| . 

For the constant barotropic model, for example, the 
explicit barotropic relation between pressure and density 



PDE^C^sPDE + {l + cl)pi,, (27) 

and the analogous k-essence Lagrangian is 0, HO] 

+ (28) 
l + cj 

where A is an arbitrary constant allowed by field redefi- 
nition in X . 



III. OBSERVATIONAL SIGNATURES 

Despite the generality of the models discussed here, 
they do provide some distinctive observational signa- 
tures. We first consider features in the homogeneous 



background properties, e.g. expansion history and equa- 
tion of state, and then in the spatial perturbation prop- 
erties. 

Consider first the jerk parameter of the expansion his- 
tory [H] 



(29) 



The importance of this parameter was first emphasized 
f73l and it was presented as one of the "statefinder" 



by m 
parameters in 



was pres 



In a flat universe, j is given by 



J 



9p, dpDE , ^ 

7:ii_D_E3 (1 + wde) ■ 

2 dpDE 



(30) 



Thus, for ACDM (including the constant pressure models 
discussed in i jllBp . we see that j = 1, independent of the 
value of r^A (52). For the barotropic models considered 
here, the requirement that dpjjE / dpoE > translates 
into a simple bound: j > 1. On the other hand, for 
quintessence models, we can write [13, 



j = 1 - -Q.de[w' - 3w(l 



(31) 



Since quintessence models generally satisfy w' > 3w{l + 
w) [37] , we see that they are often characterized by j < 1. 
Thus, accurate measurement of jerk parameter j < 1 
(or equivalently w and w') could provide perhaps the 
cleanest observational signature to distinguish barotropic 
dark energy from quintessence. One can view this as 
determining the number of internal degrees of freedom 
in the dark energy physics [13, Ull . 

Unfortunately, the jerk parameter is rather difficult to 
derive from current observational data. Assuming a con- 
stant value for j, [H derived j = 2.16^0:75. While this 
might naively seem to favor barotropic models, in fact the 
assumption of constant j strongly biases the result. Typ- 
ical barotropic and quintessence models generally have 
values of j that vary significantly with time, and cur- 
rent observations are insufficient to distinguish these two 
types of models. 

Note that as j approaches unity, it becomes difficult 
to distinguish barotropic (or quintessence) models from 
ACDM using the expansion history, since the constant- 
pressure model is degenerate with ACDM. Two possible 
observational signatures exist. In the case of barotropic 
dark energy, some of the contribution to the zero-pressure 
dark component can arise not only from dark matter but 
also from the aether component. If the dark matter par- 
ticle is detected, and its relic abundance can be calcu- 
lated from its physical properties, then one signature of 
barotropic dark energy would be an anomalously high 
observed value of Qm in relation to the dark matter (and 
baryon) calculation. 

A second possible signature is the behavior of spatial 
perturbations. Since = for the constant-pressure 
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barotropic model, one might hope to distinguish it ob- 
servationally from the corresponding quintessence model 
with = 1 on the basis of perturbation growth. (Con- 
straints on the dark energy sound speed have been ex- 
plored for constant sound speed in ^55, 56]). Interestingly, 
where barotropic and quintessence models are closest in 
their dynamics (i.e. near w' ~ ^w{l + w)), they differ 
most in sound speed, and where they are closest in sound 
speed, they differ most in their dynamics. However, as 
w ^ —1, barotropic models and ACDM become degener- 
ate, at least to linear order. We can see this by examining 
the equation for linear perturbatio n g rowth in the dark 
energy fluid in synchronous gauge [53 HI] : 

6de + (1 + w)i0DE + h/2) + ZH{cl - w)Sde = 0, (32) 

where h is the trace of the perturbation to the Friedman- 
Robertson- Walker metric, and Ode is the divergence of 
the fluid velocity. The key point is that for the case 
w = — I, there is no growing mode in Sde, since the last 
term is always positive. Thus, w = —1 models cannot be 
distinguished, regardless of the value of (this point is 
emphasized in [SS] for the case of the a = Chaplygin 
gas and ACDM). 

One can also argue qualitatively from equation (j32[) 
that the growth of density perturbations should become 
relatively insensitive to cj. in the limit where w is close to 
— 1. This conclusion is borne out by detailed comparisons 
in the w — plane between models and the observations 
[55I [5^ . These investigations show likelihood curves that 
are nearly independent of for w near —1. 

Note that barotropic models as a class can never be 
"ruled out" as long as ACDM remains viable, since as a 
limiting case constant-pressure barotropic models include 
ACDM. We have shown ways, however, in which one 
can generally constrain the allowed parameter space for 
barotropic models. That said, all observationally-allowed 
barotropic models must approach the constant-pressure 
model at high redshift, making them indistinguishable 
from ACDM in this limit (save through the theoretical 
predictions for noted above). 

One can take this further and argue that barotropic 
models "predict" a value of w near —1, in a way that 
quintessence models do not. Our argument is based on 
the upper bound on w' given by equation (jlOp . This 
equation shows that, for barotropic models with nonneg- 
ative c^, the value of w for the dark energy can never 
"loiter" at a value between w = and w = —1. While w 
can lie near for arbitrarily long times in these models, 
once it begins to decrease toward —1, equation (fTO|) puts 
a lower bound on the rate of decrease. Thus, one cannot 
have arbitrarily long periods in which w has some value 
between and —1. The opposite is true in quintessence 
models; it is easy to construct such models (trackers) 
with w roughly constant and equal to nearly any desired 
value p^ . 

In terms of our aether decomposition, the slowest rate 
of decrease for wde occurs when wae = 0. In this case, 
for example, w decreases from —0.1 to —0.9 as the scale 



factor increases by about a factor of 4 (i.e. within 1.5 
e- folds). Other choices for the equation of state func- 
tion can only produce a more rapid decrease in w. Thus, 
in barotropic models, a value of w between and — 1 
must always be a transient phenomenon, leading to the 
argument that barotropic models "predict" a value of 
w near —1 (a prediction which would have been consid- 
erably more convincing had we made it a decade ago). 
Note that a w —I attractor is also present in some 
unified models for dark matter and dark energy [sol . [60I ] , 
although such unified models are outside the scope of our 
discussion. 

Since the only viable barotropic models are those that 
scale like the dominant component at high redshift, no 
other special selection needs to be applied: in this sense 
the "bug" of not being able to distinguish a constant- 
pressure barotropic model {wae — 0) from ACDM is 
really a "feature" of ameliorating the problem of fine tun- 
ing initial conditions. (The usual cosmological constant 
problem remains of why the non- aether part, poo, is so 
small.) 

Note that the coincidence problem is also somewhat 
ameliorated (see [6l| for a recent discussion of the coin- 
cidence problem from a novel perspective). For the case 
of a cosmological constant, the ratio of dark energy den- 
sity to matter density has to increase by nine orders of 
magnitude between recombination and today, leading to 
the question of why dark energy overtakes matter basi- 
cally now. In contrast, the barotropic models can easily 
have a more natural-seeming ratio of order 10~^ or 0(1) 
at recombination. While quintessential tracking models 
can also have such ratios, they have difhculties in then 
achieving w ~ —1 today. Thus barotropic models have 
attractive characteristics with regards to both the fine 
tuning and coincidence problems. In effect, the aether 
component of the barotropic fluid anaesthetizes the cos- 
mological constant against the pain of fine tuning. 

IV. CONCLUSIONS 

Barotropic fluids have a number of characteristics that 
make them an interesting class of dark energy models. 
They have, by definition, an explicit rather than im- 
plicit equation of state relating the pressure and the en- 
ergy density. While this relation is nominally quite gen- 
eral, we show that simple physical conditions such as 
stability and causality severely restrict the allowed func- 
tional forms. In particular, we demonstrate that viable 
barotropic models must possess the following properties: 

• Asymptotic future de Sitter state, where the dy- 
namics freezes to a cosmological constant state, 

• Dynamics distinct from much of quintessence, lying 
in a separate region of w-w' phase space; sound 
speed generally distinct from quintesssence, 

• Acts as a sum of a cosmological constant and a 
perfect fluid "aether" component with wae ^ 0, 
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• The aether component must have wae ^ 0, ^ 
in the past in order not to violate matter domina- 
tion. 

These results both unify a number of special cases in the 
literature and rule out several models. 

We consider several observational signatures to distin- 
guish barotropic fluids through both the effects on back- 
ground expansion and on perturbation growth. For ex- 
ample, in the barotropic case the sound speed is not re- 
stricted to be the speed of light, as in canonical, min- 
imally coupled scalar field models. Constant-pressure 
(Cj — 0) barotropic models are however degenerate with 
ACDM, though possibly distinguishable through a dis- 
crepancy between particle physics predictions for the 
dark matter density and cosmological observations. 

Finally, the aether component of the barotropic fluid 
can anaesthetize the cosmological constant against some 



of its flne tuning and coincidence problems. In the high 
redshift universe the dark energy appears like ACDM, 
but with a dark component energy density that can be 
comparable to the matter density. At late times, it nat- 
urally and rapidly transitions from a matter-like behav- 
ior to behavior that approaches a pure cosmological con- 
stant. 
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